Negative differential conductance induced by spin-charge separation. 
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Spin-charge states of correlated electrons in a one-dimensional quantum dot attached to inter- 
acting leads are studied in the non-linear transport regime. With non-symmetric tunnel barriers, 
regions of negative differential conductance induced by spin-charge separation are found. They are 
due to a correlation-induced trapping of higher-spin states without magnetic field, and associated 
with a strong increase in the fluctuations of the electron spin. 
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In single electron transport 0, the electron charge 
causes the Coulomb blockade effect that stabilizes the 
electron number in a device. Experiments show also that 
the electron spin can lead to important mesoscopic trans- 
port effects. For instance, quantum dots in carbon nan- 
otubes [3-0- 01 have revealed several non-equilibrium and 
coherent spin processes, and spin-parity effects. 

For few-electron quantum dots, the spin blockade effect 
has been predicted p|- In one dimension (ID) this is indi- 
cated by a negative differential conductance that occurs 
only when a state with maximum spin value is occupied. 
Combining spin blockade with spin-polarized detection, 
the electron spin in a 2D-lateral quantum dot was probed 
Also, in a few-electron vertical quantum dot higher- 
spin states, not directly accessible from the ground state, 
were detected via probing with non-equilibrium voltage 
pulses 0- When spin relaxation was absent, strong fluc- 
tuations of spin and charge were observed. Motivated 
by experiments on a quasi-ID quantum dot formed by 
two impurities in a quantum wire Q, we have calcu- 
lated microscopically single electron transport through 
a ID quantum dot formed by two equal tunnel barri- 
ers in a Tomonaga-Luttinger liquid. Signatures of non- 
Fermi liquid behavior and, in the presence of a magnetic 
field, correlation-enhanced spin polarization were found 
@. Spin-charge separation has been observed and 
analyzed theoretically [TfT | in tunneling between parallel 
quantum wires. In the coherent tunneling regime, the 
effects of non-Fermi liquid correlations in ID quantum 
dots have been investigated non-perturbatively |l2|, |lj| • 

In the present work, we consider a ID quantum dot 
in the incoherent sequential tunneling region. We pre- 
dict that states with higher total electron spin in a non- 
magnetic ID quantum dot containing an arbitrary num- 
ber of electrons can be dynamically stabilized by tuning 
the asymmetry of tunnel barriers. The occupation of 
these states can lead, in the sequential tunneling regime, 
to a negative differential conductance. This is caused 
by the peculiar nature of the non-Fermi liquid correla- 
tions that determine the tunneling rates. The new phe- 



nomenon is due to spin-charge separation and is different 
from the previously discussed negative conductance asso- 
ciated with spin selection rules via Clebsch-Gordan coef- 
ficients that were included in the tunneling rates ad hoc 
0- We study the combined effects of asymmetry, elec- 
tron correlations and relaxation. We predict that sta- 
bilization of the higher-spin states is associated with a 
strong increase of the spin fluctuations. Relaxation of 
these states strongly affects the negative differential con- 
ductance, in contrast to "normal" states associated with 
positive conductances. Our results open a novel possi- 
bility of experimentally addressing non-Fermi liquid be- 
havior without performing often cumbersome analyses of 
temperature dependences. 

We start from the microscopic Hamiltonian H = + 
H\ + H t , which contains the dot (H^), the leads (H\), and 
a tunneling term (H t ). The dot is modeled as an inter- 
acting ID system confined to \x\ < d/2, with excitations 
associated with independent energy scales 



H d {n,s,p,a) = -j-(n-n g ) + —s 
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Here, n and s represent the number of electrons, mea- 
sured relative to the number of charges n g corresponding 
to a gate voltage V g , and the z-component of the total 
spin in units of ?i/2, respectively. Energy scales E p and 
E a are the charge and spin addition energies. The last 
two terms do not change the total charge and spin but 
describe intra-dot charge and spin density waves, where 
p and a are related to bosonic creation and annihilation 
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operators via p = J^j 3% j^p-j an( l G 

The energy parameters in Eq. reflect the micro- 
scopic electron interaction parameterized by |l4j 
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with the exchange and Coulomb matrix elements V x = 
V (2kp) J '2ttHvf and Vq = V(q = 0)/2ttHvf, respectively. 
In experiment, the charging energy is not solely influ- 
enced by the microscopic interaction in the quantum 
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dot. Therefore, we treat it as an independent param- 
eter, typically E p 3> E a . The spin addition energy, 
E a = ■nhv a /2dg (J1 is due to the Pauli principle. The ex- 
citation energies of the neutral charge and spin modes 
are e v = Trhv u /d with charge and spin mode veloci- 
ties v u = «f(1 + V^)/g v [y = p,a). Without inter- 
action, g p = g a = 1, we have E a = Eq = TrhvF/2d 
and s p = e a = 2Eq. We assume the leads to be Lut- 
tinger liquids with Coulomb repulsion go < 1, low-energy 
charge and spin mode dispersions oj p (q) = v-p\q\/go and 
ui a (q) = vp\q\. We assume go ^ g p , since the interaction 
can differ between dot and leads. Tunneling between the 
leads and the dot at x = ±d/2 is described by H t with 
amplitudes Ai and A r . We consider weak tunneling 9]. 

For not too low temperatures, k B T 3> SE, transport 
is dominated by incoherent sequential tunneling and one 
can safely neglect higher order coherent processes [lij . 
Here, SE represents the level broadening due to higher 
order contribution from tunneling via virtual states, and 
is proportional to the rate 7 in Eq. The kinetic vari- 
ables are then n and s only since p and a relax towards 
thermal equilibrium very efficiently due to mechanisms 
such as coupling to phonons and spin-orbit interaction. 
The opposite limit, but for the spinless case, has been 
considered in The relaxation of the spin s, which is 
effected by the same processes as well as co-tunneling, is 
usually much slower |l7]|. Therefore, both tunneling and 
spin relaxation have to be included in the Master equa- 
tion. In order to be able to resolve the spin dynamics we 
consider k-g,T < E a . In the stationary limit, one has to 
solve J2 V > [Piv')^'^ ~ P(v)r n ^ r '} = for the probabil- 
ities P(rf) = P{n 1 s). Tunneling is characterized by rates 

r^'' , subject to the selection rules n' = nil, s' = s±l. 
These contain the non-Fermi liquid correlations and can 
be calculated microscopically for finite temperature @. 
For simplicity, we quote only the expression for T = 
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with uj c the cutoff frequency ^3 1 Gr,i = (7reA r j/[j c ) 2 /h 
the tunneling conductances, and E-^ obtained from the 
charge and spin addition energies corresponding to the 
states i and j. The neutral spin and charge modes (en- 
ergy scales e p and e CT ) are fully taken into account. For 
the spin relaxation, obeying the selection rules s' = s ± 2 
and n' = n, we use the detailed-balance result 
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FIG. 1: Differential conductance at ksT — 10~ 2 E IT in the 
(V ,n s ) -plane with asymmetry A = 100, interaction strength 
in the leads go = 0.9, charging energy E p = 25E a , and spin 
relaxation w = (V in units of E a /e)\ (a) e p = 3.0-E CT , s a = 
2.2E a \ (b) e p = 15-E CT , e a — 2.5Ecr; (c) states involved in the 
transitions for V > near the linear conductance peak, and 
denoted by the triples of integer quantum numbers (s,a,p); 
(d) non-interacting case, e p = E a = 2£"o. Right edge: grey 
scale in arbitrary units, labels and scales in a,b,d are identical. 



containing a phenomenological rate u>, and 1//3 = k^T. 

Numerical results for the differential conductance G in 
the plane of bias V and gate voltage V g oc n g are shown 
in Fig. [2 for g = 0.9, w = 0, k B T = 10~ 2 E a , and non- 
symmetric barriers near a charge transition n <-> n + 1 (n 
even). The asymmetry A = G\/G r has been assumed to 
be large, and the temperature low, in order to emphasize 
also weaker features in the conductance. The intersection 
points at the centers of the grey-scale panels (n g = n ICS , 
V = 0) denote a linear conductance peak. The grey 
areas denote the Coulomb blockade regions. Black and 
white lines for V ^ correspond to positive and negative 
differential conductance peaks, respectively. Here, tran- 
sitions involving additional spin states as well as neutral 
spin or charge modes in the quantum dot become relevant 
(Fig. ^) . Remarkably, it is the spin-charge separation in 
the dot that induces negative conductance peaks which 
separate certain spin states. Independently of the asym- 
metry A, these features are not present when spin-charge 
separation is removed, i.e. e p = e a = 2Eq (Fig.^Ji). 

In order to understand the physics behind the negative 
differential conductances we concentrate in the following 
on the regions denoted as I, II and III (Fig. Here, 
for k^T < E a , only the five states 77 = (n,s = 0), (n + 
l,s = ±1), (n, s = ±2) are necessary to obtain a good 
approximation for the conductance. We assume A > 1 
and that the electrons flow from right to left for V > 0. 
Then, the quantum dot with lower particle number n 
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will have a higher occupation probability. This implies 
that negative conductances can only occur parallel to the 
transition line (100) -> (000), see Fig-Ufc. 

Using this model, the Master equation can be solved 
analytically using the rates from Eqs. ©-© and the en- 
ergies Ea-,i = eV/2 + E p (n g -n l . cs ), E 2 ^i = E ^i + 2E„ 
and Ei^q = eV — £?o->i) Ei_ 2 = eV — E 2 ^i- The cur- 
rent /(V) can be evaluated, and from that G = 81/ dV = 
N/D. Keeping only the dominant terms near the transi- 
tion lines parallel to (100) — > (200), with go « 1, one ob- 
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From Eqs. (JBJ) and J7J one recognizes that for obtaining 
a negative differential conductance one needs N < 0. 
This implies 8J2 p K { 2 p) < -^T, P K [ P) since tfEpA^ 
is always positive. From Eq. |j3J one finds that the rate 
changes depend crucially on the presence and strengths 
of the correlations via a v . Without correlations (go = 
9p = H<t = 1), the rates are integer multiples of 7 r j. In 
particular, one finds T 2 ^ 1 = 2T® such that always 
G > (Fig.nji), even in the most favorable limit A — > oo 
and w = 0. This suggests that it is the influence of the 
intra-dot non- Fermi liquid correlations on the rates which 
yields the negative differential conductance. 

To support this we introduced correlations in the dot 
(g a > 1, g p < 1) while keeping the leads non- interacting 
(go = 1), for w = and A — > oo, with a fixed G r . One 
finds that along the line (100) — > (200) one has A^ 1} > 0, 
and Aj ^ = 0. This implies always negative conduc- 
tance in I and II since 8 < 0. The state (200) acts as 
a bottleneck for the current, since it cannot relax to the 
ground state via a single tunneling event. It accumulates 
probability at the expense of (000). In region III, due to 
the additional activation of the charge wave state (101) 
in r 2- * 1 , depending on a CT=1 + a p=1 1 one has 8^1 
and the conductance can be both negative and positive. 
However, for interactions with Vb > 2V X one always gets 
G > 0. In Fig. \I]p, region III is absent due to g p <C 1 
(strong correlation). Therefore, we find only negative 
conductance as long as (300) does not contribute. Sim- 
ilarly, one can analyze the conductance associated with 
higher quantum numbers (s,a,p). In all of the cases ana- 
lyzed we have found that the necessary ingredient for neg- 
ative conductance is the spin-charge separation. In order 
to show that the states with higher total spins are re- 
sponsible for the occurrence of negative conductances, we 
have studied the influence of spin relaxation with vary- 
ing strength. Figure [21 shows specific numerical results 




FIG. 2: Differential conductance (units e ■y T /E tT ) at positive 
voltages V, in units E a /e (right), near (100) «-» (200), and 
at —V (left), for g p = 0.135, g a = 1.25 and go = 0.9, at n g 
in region I, with k^,T = 10~ 2 -E CT . Top panels: w = 0, with 
A = 1 (red), 2 (green), 2.6 (blue), 5 (magenta) 50 (cyan). 
Middle: A = 10 with w/j = 0.3 (red), 0.35 (green), 0.4 
(blue), 0.45 (magenta), 0.5 (cyan). Bottom: phase diagram 
of the crossover between positive and negative conductance at 
T = for g = 1 (red), 0.9 (green), 0.8 (blue), 0.7 (magenta), 
0.65 (cyan) for n g in region I (w in units 7 = r )\(E a /huj c ) a ). 



for the differential conductance at positive and negative 
biases for n g fixed in region I near the (100) — > (200) 
transition, when the asymmetry A (top panels) and the 
relaxation w (middle panels) are varied. 

At w = 0, for small asymmetries, say A < 2, all peaks 
are positive. For larger asymmetries, A > 2, the peak 
corresponding to positive bias can become negative (top 
panels). At fixed asymmetry (A = 10, middle panels), 
by increasing the relaxation rate of the state (200) to- 
wards the intrinsic rate 7 = 7i(A (T / fiw c ) a , the height of 
the negative conductance peak is reduced, while the cor- 
responding positive conductance peak for V < is un- 
changed. When w -C 7, the (200) state is stable and the 
negative conductance feature is strong. However, when 
w > 7/2, the (200) state becomes unstable in favor of 
(000) and the negative conductance vanishes since the 
ground state channels are re-opened. 

For finite asymmetry, A > 1 , and with spin relaxation, 
w > 0, one finds a finite critical value, A c , for each neg- 
ative conductance feature such that one gets G > for 
A < A c . The critical trajectories w c (A) for which the 
negative conductance peak disappears depends on the 
state with the highest spin s max > 1 involved in the trans- 
port, on the interaction strength in the leads, and on the 
relaxation. In the above example with s max = 2, w c (A) is 
obtained by setting N = in Eq. JJj} at the voltage cor- 
responding to the position of the negative conductance 
peak and finite interaction in the leads. This phase tra- 
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FIG. 3: Differential conductance (red curve, left scale, units 
e 2 7 r / ' Ea) and spin variance (blue curve, right scale) as a func- 
tion of bias voltage V (unit E a /e) for n g in region I, A — 100, 
w = 0, k B T = 210~ 2 E a , E p = 25£ CT , e a = 2.5E a , e p = 15E a , 
go = 0.9. 

jectory is shown in Fig. (bottom) for different go- The 
negative conductance depends crucially on the interac- 
tion in the leads go ■ For decreasing go tunneling becomes 
less efficient, giving spin relaxation increased importance. 
However, with interactions in the leads, the negative con- 
ductance features, once they arc stabilized, can even be- 
come enhanced as compared with their positive counter- 
parts at negative voltages. Similar results are obtained 
for higher s max , but spin relaxation is then more efficient 
in preventing their formation, w c (s max + 1) < w c (s max ). 

Figure [3] shows the differential conductance and the 
spin fluctuations 8s = E n s( s — (s)) 2 P(n, s)] 1 / 2 as a 
function of the bias for go = 0.9 with A = 100. The con- 
ductance shows positive and negative peaks for increas- 
ing bias (beyond the first two ground state peaks not 
shown here) . One observes a step in the spin fluctuation 
near the negative conductance peaks. This is due to the 
participation of the higher-spin states in the transport. 
For example, the first negative peak at V = 2.3E a /e in 
Fig. El contains mainly contributions from s = 2. Con- 
sequently, the spin variance jumps to a value dominated 
by P{n, s = 2). This remains dominant even if new con- 
ductance channels enter (subsequent positive peaks). It 
changes only when the s = 4 spin channel enters (second 
negative peak at about V = 7.3E a /e). Since we consider 
the transition n <-> n+ 1 with n = even, we conclude that 
for positive voltage, even values of the spin dominate the 
spin dynamics. For negative voltages odd spins are dom- 
inant. Strong steps in Ss are also found at the voltages 
of the positive conductance peaks with n g in region III. 
This indicates that they are an intrinsic feature of the 
contribution to transport of the states with higher spins. 

In conclusion, we have investigated the non-linear, se- 
quential transport in a ID non-Fermi liquid quantum dot 
embedded in a correlated electron system including spin. 
We have found that there are distinct negative differen- 
tial conductances induced by spin-charge separation, if 
asymmetric tunnel barriers connect the quantum dot to 
the leads. We found that one can - without applying a 
magnetic field - stabilize states with higher total spins in 
certain regions of the parameter plane spanned by bias 



and gate voltage if the asymmetry of the tunnel barri- 
ers exceeds a critical value depending on spin relaxation 
mechanisms competing with the non-Fermi liquid corre- 
lations. The participation in the transport of the states 
with higher total spins is indicated by the sensitivity to 
spin-flip relaxation and by strong changes in the spin 
fluctuations. The predicted phenomenon should occur in 
quasi-lD quantum dots containing even many electrons 
for spin states different from the ground state but not 
necessarily with maximum total spin. The physical ori- 
gin of this new correlation-induced trapping phenomenon 
are the non-Fermi liquid properties of the system which 
lead to spin-charge separation and non-trivially enter the 
tunneling rates. The results in Fig. show that the neg- 
ative differential conductances are not restricted to very 
large asymmetry. Depending on the interaction param- 
eters, they can occur also for moderate asymmetries of 
the order of 2 and even smaller. We expect that the 
above non-Fermi liquid phenomenon will be accessible 
using quasi-lD electron systems with moderately con- 
trollable tunnel contacts such as carbon nanotubes, and 
quasi-lD semiconductor quantum wires. For the latter, 
asymmetric barriers seem to us to be the genuine case. 
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